et Popular Computing 34 


January 1976 Volume 4 Number 1 


A Problem for Backtracking 


PC34-2 


.5314789170422551237539087 890528 300567 75757259887155 
526 36052461616 13896667667 39 33130310366 3703146946000 
-83095189484 5300470874 152877 54 558 3076521 3983 34885972 
. 23961180127 748338407 146992427 20297 00378378968526529 
.02439745849988504.25108172455419374.1911462170107 3118 
.42281321982187 2839528117 387131165958496133027248787 
0358927 39572482224 3970172784810137947509777064 39141 


e34 583461742527454 .8814029027 346103910190036592389411081 
057829421204 31667674211950581147 10386 


# H inv) Ww WwW Ww bed 
NESERIES. 34 


Aah 80001047 150456 339 .5529492519141529488801881667984087 3 
60697452457 315 3636494900997 26723166 


tan7+ 34 1.541393038590891503957 2188274 16327 28484 3325036165643 @ 


PRINTED BY 
STATI 
Pp) panier Np, 
‘argeld 


17002 Ventura Blvd, Encina, Ca. 91316 


POPULAR COMPUTING is published monthly at Box 272, Calabasas, California 91302. Subscription 
rate in the United States is $18 per year, or $15 if remittance accompanies the order. For Canada and 
Mexico, add $4 per year to the above rates. For all other countries, add $6 per year to the above rates. 
Back issues $2 each. Copyright 1976 by POPULAR COMPUTING. 


Publisher: Fred Gruenberger Contributing editors: Richard Andree Advertising Manager: Ken W. Sims 
Editor: Audrey Gruenberger William C. McGee Art Director: John G. Scott 
Associate Editors: David Babcock Business Manager: Ben Moore 


Irwin Greenwald 


ey Reproduction by any means is prohibited by law and is unfair to other subscribers. ee 


@ 2023 This work is licensed under CC BY-NC-SA 4.0 


PC34-3 


@ Backtracking {| BY THOMAS R PARKIN 


It is suggested that the serious student of this backtracking 
topic provide himself with a set of cut-out squares corresponding 
to the markings on the problem on the cover, in order to aid in 
following the discussion. 


Last month we saw that backtracking is a simple 
procedure, in concept, and often quite difficult to apply, 
largely because of the two requirements to generate the 
solution space of the problem systematically and to be 
able to detect the undesirability of further exploration 
along some branch at an early point. We shall apply 
backtracking to the problem on the cover as an illustration 
of how it works and where the difficulties lie in its 
application. 


The problem consists of an arrangement of 16 squares, 
each capable of being placed in any one of four rotational 
orientations in any one of 16 locations in the larger 

2 square. We thus have 16! ways to locate the 16 small 
squares in the larger square, and 4 to the 16th power 
ways the individual squares can appear due to rotations 
for each way they are located. This gives us 


y26 361 or about 3.5°1018 


possible arrangements comprising the solution space. The 
end case (or cases) of interest must satisfy the condition 
that at each of the interior points of the larger square, 
where 4 small squares meet, such points must be surrounded 
by 4 identical numbers on the corners of the small squares. 
Now, even if we could construct the individual cases and 
test them in one microsecond each, it would still take about 
100,000 years to try every case! So, let us apply some 
analysis and backtracking to see if we can cut that time 
down. 


We shall leave internal coding details until later 
and concentrate on how we should go about exploring the 
solution space in a way that will allow us to throw out 
very large numbers of cases without considering them. 


First, let us note that there are seven of each 
corner number from 1 to 9 except for the number 4, of 
@ which there are 8 each. Next we note that there are 9 
places where four corners meet, so we either have all the 
numbers at one each of the corners or all but one at seven 
corners and the number 4 at two corners. 
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Now, we could start by picking any one of the 16 
squares, rotate it to any one of its 4 positions, and 
place it down in the, say, upper left corner. Next we 
would search among our remaining 15 squares for one 
(or Bones which would match that square at the corresponding 
corners (lower right for the first square; lower left for 
the second square). We would find one or more such 
matching squares at this stage of the process, of course, 
but maybe we won't be so lucky later. We now continue 
this process (i.e., in an orderly way), progressing from 
left to right, top row first and then successive rows, 
selecting from our remaining store of squares one which 
matches at the corresponding corners according to the 
conditions of the problem. Place that square in position 
and continue. Thus we have our systematic generation 
mechanism. 


Sooner or later we shall either come to a halt 
because there is no remaining square which can be selected 
to match in the next place where we have a square to go 
(this is our test mechanism} or we shall use up all 16 
squares and we have a solution to the problem. 


Let us look at this process in more detail. Suppose 
we simply start with square A where it is given in the 
diagram and look for a square to match it on the right. 

We can now pick from ©, F, H, L, FP, or Q, each of which 
has a 6 to match Ats 6 in the lower right corner. Let's 
just take this next square in alphabetical order (which 
is the rule we will follow all the time), or square C. 
The process continues this way, then, until we have 

A, C, E, B across the top row, with C and E each rotated 
90° counterclockwise. F now goes in the first position 
of the second row. At this point, the test for a matching 
square gets somewhat harder--there are three corners to 
match simultaneously: the 6's corner, the 7's corner, and 
the 8 corner (on F). We have the situation shown in 
Figure l. 


Figure 1 


We start looking through the remaining squares for one 
which will match at the 8, 6, and 7 corners. The only 
candidates (with at least the 6) are H, L, P, and Q. H 
comes close, since at least both the 8 and the 6 match, 
but not the 7 corner; there, H has a 1, and hence won't 
do. None of L, P, or Q will work either. Well, here 
we are at a road biock. We still have 11 squares to use 
(representing 11!-4 to the llth power possible cases) but 
we know that not one of these can possibly be a solution if 
the first 5 squares we have placed in the first 5 locations 
as shown remain exactly in that arrangement. We have thus 
eliminated a very large part of our solution space simply 
by picking five squares in some logical systematic order 
(in this case, alphabetically), in the order of potentially 
eligible squares for each next location. 


Now we backtrack. Let us undo the last step we 
took prior to reaching the impasse; that is, backtrack 
one level, discard the choice we made (F), and choose the 
next eligible square for that location and try again to 
continue. This results in trying H in the first location 
of the second row. Again we are blocked, because among 
F, L, P, and Q none fits properly to match at the 3 
corners where it must match. Again, we backtrack one 
level and try L. Again we are blocked. We continue 
this process of backtracking, proceeding until blocked, 
and repeating the process until we have exhausted all the 
squares that will possibly fit in this fifth spot. Of 
course, if one fits and we can put in still another that 
fits, we simply continue the process until we either have 
a solution or are again blocked, whereupon we repeat the 
back up-try again cycle. In this case, we shall find 
that none of F, H, L, P, or Q can be used as the first 
square of the second row (always holding the first row 
constant as we have selected it) and allow us to pick any 
square to fit in the second spot in the second row. 


We are now in the position of having exhausted all 
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possibilities for the 5th square, so let's backtrack again 


to the 4th square. 


At this point we have now proven (note that this 
process constitutes an actual proof) that with the first 
four squares fixed as shown in Figure 1, there is no 
possible solution which satisfies all the conditions of 
the original problem. Let us now try the next eligible 
square for the 4th position. Possibilities are I, J, L, 
and @ and, of course, we pick I, since we always select 
from among those eligible in alphabetic order. Then, we 

o to the 5th position where the candidates are, again 
but not necessarily) F, H, L, P, and Q. At this point 
we have a choice--either continue the process exactly as 
we have been doing (in which case we will find that no 
combination will allow us to pick both a 5th and 6th 
square) or we can get just a bit clever and, having 
remembered what we did, simply reject the choice of I 


b 
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for the 4th place and backtrack again immediately. In 
order to take this latter step, we have introduced the  ) 
complication in our program of storing previous sequences 

of actions. In a simple case like this one, perhaps that 
refinement would not be justified, but the idea comes in 

the category of tricks to shorten running time which could 
be very valuable in some more extensive problem. 


In this sample problem, we will discover that we 
run through all of the possible 4th squares without ever 
being able to get a 6th one in place. So, again we 
backtrack to the 3rd square, and repeat the continuing 
constructive process. (At this point it is no fair 
indicating what will happen-~elther try it by hand or 
write a computer program to solve the problem.) ssi 
perchance, we continue this going forward and backtrack- 
ing process until we get to the point where we must change 
the first square, we must now remember that we simply 
placed it down in the orientation it happens to have in 
the original statement of the problem on the cover. 
Hence, we must rotate that square by 90°, and proceed 
as before. If we should find ourselves back at level #1 
again, again we must rotate, and, of course, still again, 
until we have tried the first square picked in all its 
possible orientations. Indeed, if at some point we have @ 
found more than one matching number on, say, the 3rd 
(or any other) square, we must remember to try the second 
possible orientation of that square before we discard it 
in that spot. 


However, suppose we have exhausted square A in all 
its four orientations in the first place, then we must 
start the process anew with B and so on. If we follow 
this procedure through all 16 squares, each in all four 
orientations, we will have totally exhausted the solution 
space of the original problem without having actually 
looked at more than an extremely small fraction of the 


possible end cases. Furthermore, in the process we 
will have found all of the solutions (if any) to the 
original problem. We could also have stopped after 


finding the first solution (if any) if we wished. 

Frequently in practical combinatorial problems, ali, or 

at least the number of, solutions are wanted, and if 

there is none, running through all cases is necessary to 
demonstrate that there is no solution. The application 

of backtracking through the use of suitable tests to a 
suitable orderly generation of cases can very materially 
shorten the solution time for a great many problems. 

Indeed, even in this simple case, a program running a few r ) 
seconds to a few minutes can solve this problem which » 
would be unthinkable to do by brute force. 


While it might not be necessary, for this particular 
problem, to be terribly clever about the lexicographical 
order of generating the solution space, perhaps a further 
point could be made to illustrate that this process is 
subject to refinement and different approaches, which can 
be powerful aids in shortening the running time for some 
problems. 


For example, it is obvious that the requirements of 
this problem place the most restrictions on the choices 
of the four center squares. These must not only match 
in the center, but at two other places for each square. 
Thus, if one were to start the solution by picking one of 
the four center squares and then trying to fit 3 more 
around that center point, it turns out to be much tougher 
to get even the 3rd or 4th square into the pattern, and 
the blockage at any point excludes even more cases at 
each application of backtracking. However, the coding 
and programming become quite a bit more difficult, and 
perhaps too hard to implement for their advantageous use, 


We are now in a position to produce a flowchart 
for this problem. Figure 2 is the flowchart at a fairly 
gross level showing the basic logic we have used. It 
is at this point that coding details enter the picture. 
For example, how should we represent the squares for 
internal handling by our program? How do we remember 
where We are in each list of potential next squares for 
each position? To what extent can we trade storage space 
for running time? Should we maintain lists or recreate 
them each time? Answers to these and a host of other 
questions will determine how we design our detailed flow- 
chart and how we program the problem solution. For 
example, it becomes quite clear shortly that each square 
can be thought of as a set of 4 triples; that is, the 
square consists of four corners, each having a digital 
value, and each corner has a left and a right corner 
corresponding to it. We will find that we are concerned 
about the two values on either side of a corner of a 
square when it is time to match that square into position 
as the next square in sequence. Hence, perhaps we should 
adopt an internal code representation of a square which 
looks something like Figure 3. 
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such as (A, 6, 9, 4, 8), but remember that we must be 
able to remember the orientation of the square when we 
put it down and it must match in from 1 to 4 places with 


other squares. 


squares a bit easier to work with. 


We may find the redundant coding of the 


Similar questions of coding and indexing of work 
space and provision for lists and counters and pointers 
will comprise the preliminary considerations for our 


program. 


ease of creating the program, 


While these details can materially affect the 
they are merely details. 


The essential items to be thought through when potentially 


applying backtracking to a problem are these: 


1) Devise 


an orderly procedure which will generate the totality of 


end cases in the solution space of the problem, 


and 2) 


Devise the most general criteria or tests which can detect 
dead-ends in the generation process at the earliest 


possible time. 


With these two tough (usually) questions 


answered, the application of the conceptually simple 
process of backtracking becomes easy. 
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The CSR Function e 


aq,+ ¥o+ ya,4 da, + 4a + ¥ a veneme 


The function A, named CSR for "continued square root," 
is most interesting. 


If the successive a's are: 


eo wen 20 mts bbe 1... 89a. LOO sna eee 
(successive differences are 4, 6, 8, 10,...) 


the CSR converges to precisely 2. 
But if the a's are the factoriais: 


ly Be, B BB, 120, 720, 5040, 40320, 362880,... 
(which are vastly larger) 


the CSR converges to the following (computed by Dorothy Cady): 


1.82701 47176 08592 22637 38431 92852 89247 37479 36296 
re ane 61624 62956 21001 52387 09672 78310 72441 
3 055 


This is fun. If the a's are all one, the series converges 
to the Golden Mean, (1 +/5)/2. 


Other values to try for the a's are the following: 


Consecutive integers starting with 1. 
Consecutive odd integers starting with 1. 
Consecutive even integers starting with 2. 
Fibonacci numbers: 1, 1, 2, 3, 5, 8, 13,... 
Squares, cubes, and higher powers. 

Prime numbers, starting with 2. 


WN 
= 
= 
= 
a 
2 


amo Q0 > 


> 
@ 
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Herman P. Robinson has shown that the a's can be 


chosen to converge to any positive number, if both 
negative and positive terms are used, and that there are 
infinitely many ways in which this can be done for any 
number desired. To produce the integer A, one method 


could be: 
a, =i 


an 


with K arbitrary 


(orn = 2K 0a" @ (nee ee ee 


For A = 2, K = 1, we have the sequence given above (at #) 
that converges to exactly e. Ione fh 


sequence of a's is the following: 


Dy We By wae 


lp dig Sikh Ui, soil, Sol, 7el, MOS, BON, UsOl,... 


The sequence of a's can be made to converge to 
irrational (or transcendental) numbers, 
various constraints applied to its values. For example, 
the following sequence contains only prime numbers and 


converges to7z7: 


e 5 17 37 5S 
317 859 887 1637 
6491 10532 13399 14083 
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Words of storage from address OO0O1 up are cleared 
to zero. Starting at word 0001, each word of storage 
is to be loaded with a number, according to this scheme: 


If the contents of a word is non-zero, leave it alone. 


If the contents of a word is zero, it is to be 
supplied with a number. If its address is A, then link 
it to the word at address B, where B = 3A+1. Word A 
is to be supplied with the number B, diminished by the 
amount in the first word to the left of word A. Word B 
is to be supplied with the number B, increased by the 
amount in the first non-zero word to its left. 


The contents of storage will appear as follows, for 
the first few words: 


4 3 i 8 8 (o) 15 0 0 _25 
ese fk  5 6, Te me 10 
9) @) 0) 0 ) 41 Oey pm On <5) Oe 
oa “2 "43 ay aG ae 1" 18 19 20 


At the stage shown, words 1 through 5 (and their associates, 
words 7, 10, and 16) have been supplied. Word 6 is next, 

so at this point A = 6 and B = 19. Word 6 is supplied with 
19 minus 8, or ll. Word 19 is supplied with 19 + 41 = 60. 


We now have some problems: 
1. What will be the contents at address K? 


2. Given the contents of a word as J, is its address 
uniquely determined? 


prosiem 113 


3. After the process has passed L words, will any 
word at an address less than L have zero for its contents? 


4, If the maximum absolute value that a word can hold 
is M, how many words can be filled before the process must 
terminate due to overflow?» 


O 
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Send us $100 per month 
for 12 years and we'll 
send you $100 per 

month forever 


A financial institution advertises: 


This is simply a clever way of expressing 
confidence that the institution involved 
can maintain a 6% interest rate 
(compounded monthly) indefinitely. At 
that interest rate, $100 per month for 
144 months builds up a reserve of over $21000, and $21000 
at 6% can then generate $100 per month in interest without 
disturbing the principal. As a matter of fact, the 
build-up period could be 139 months, rather than 144, and 
the same offer could be made (the amount available at the 
end of 139 months would be $20004.84), 


They could also change the offer 
to read: ——-———_——— 


Send us $100 per month for 
l2 years and we'll send you 


20 th 
e Again, the same offer could be made using ee pee oa tee 
139 months instead of 144, 


All the above assumes constant interest 
rates, and the problems could be solved quickly 
and readily with compound interest tables. 


Suppose, however, that the interest rate were not 
constant. Assume that a fund of exactly $20000 is 
available for payout, at an initial rate of 6%, compounded 
monthly, but that the interest rate rises by .00008333333 
per month; that is, the rate rises to 7% at the end of ten 
years and continues to rise at that rate. How much can 
be paid back per month, so that the fund is completely 
depleted after 139 months? 


PROBLEM 1 14 


PC34-14 


Fibonacci 


In issue 25 and again in issue 30, selected pairs 
of successive terms in the Fibonacci sequence (numbered 
so that the 12th term is 144) were given in full. The 
following list gives successive pairs to 50 significant 
digits, plus the proper exponent. 


1500 1. 355112566856 31019510 369 36867 148408377786010712418 
1501 2.1926181917556241406686103706 309915958486962 357678 
2000 «4.224696 3333923048787067256023414827825798528402507 
2001 6.835702259575806647045 396549170580107055408029 3655 
2500 1.31709051675194 962952276 30871253164 120666069649925 
2501 2.13109722236481725896 324299517047974 31818205 363701 
3000 §=4.10615886307971260333568 378719267 105220125 108637 37 
3001 6.64390460366960072280217847866028 38424416 351245278 
3500 = =1.280135297794681 361535851 3682510196153890048112207 
3501 2.071302422030262753416199537 299 34 1464207691967 3446 
4000 3.990947 34 35004422792081248094960912600792570982820 
4001 6.457488449094817 353137694901536959564441 3900640151 
4500 1.244216976598660009 311270917708920 361214744 3446630 
4501 =. 2.01. 3185 3575162644 311220549 3352144 34180594 361169251 
5000 3.878968454 3883256 33701916 3083259053120821277 146462 
5001 6.2763028004889570860 35253108 3496840554785287027 365 
5500 =: 11.209 30645965564.4637 19333621498229900647 52902309925 
5501 1.9566989545376364830250428708616122226106623318352 
6000 =3..77013149387799945333900224880775791236117 33793010 
6001 6.1002008991510192252609003226598130375096219 309442 
7000 3.664348305037232832276358967281604921857 154 3934176 
7001 5 .929040104168 3102440672829867913655667 177284295571 
8000 = 3..5615332004606267 397689149054 27460387 1413695391102 
8001 5.76268177687862267446233299765 388917247 12626724313 
9000 =. 3. 46160291 286684746 31328927 294.06531958210049 3884057 
9001 5.6009911685741999046494 37601317 3821341097 357445056 


10000 
10001 


11000 
11001 


12000 
12001 


13000 
13001 


14000 
14001 


15000 
15001 


16000 
16001 


17000 
17001 


18000 
18001 


19000 
19001 


20000 
20001 


Fm FM FH FM FN FN FO FW WW WwW WW 


. 36447648764 3178 3266621612005107 54 3310 302148460680 
-443837 31135652813387 3426099 3750380135 389184554696 


+27007 5257 2827513967 58178170577 984 303744 5034881959 
-291092912053548874786829 39 367 34796 3292007 24429812 


. 1783227 576613539 368 3233 3682339347 21484745251 12428 
.14263424911336592579 396579289954 5208268 34.44 352683 


.0891446701021944 14086967 697539 39461777479 35397213 
- 9983410723909 3166884475 32760487 3161500 35287 26517 3 


.0024687611784610909954941797150256486927479 374908 
-8580965057465408 3466920629747 2292420409 1728616391 


-918224824 20491 3830236407 2236985132022309626557118 
-721786952377237415507 76991928244 59449 357584 285914 


. 83634462237 111784 31 301746147699818279 332301266922 
. 589 3020028044.540397867604 170977 3581819153596 36710 


.7567618334702578264 3154506617507 362025276807 78087 
46053434544 33546469 35622154 307 3646354422954999890 


-6794119961787633376993101361978426649041765169298 
~33537967968144 24 56 3578282567 3254 38165555868749 331 


60423245784 32 3240223747 64153109917871045715829708 
«2137 36631396027 7062669 329 31928527 577 3791917287840 


. 5311623237 3236124.2240155003520607 291766 3564858025 
.09550667084.2125091974920 36759 33507499284074 326584 


Further entries in this table are solicited, as 


follows: 


1. Even 1000's. 


E2089 
E2089 


E2298 
E2298 


E2507 
E2507 


E2716 
E2716 


E2925 
E2925 


£3134 
E3134 


E3343 
E3343 


E3552 
E3552 


£3761 
E3761 


E3970 
E3970 


E4179 
E4179 


2. Pairs of successive entries to allow for 


interpolation. 
3. At least 50 significant digits. 


PROBLEM 1 1 6 
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Problem 47 (PC14-5) was the following: 


Four blocks of storage are addressed at A through 
A+9; B through B+9; C through C+9; and D through D+29, 


Blocks A, B, and C contain numbers which are in ascending @ 
order within each block; there are no duplicates among 

these 30 numbers, We want to merge the 30 numbers into 

block D. The main routine has already verified that 

the numbers in blocks A, B, and C are as stated. Draw 

a flowchart of the logic of a subroutine for the merge. 


PC34-16 
Problem Solution 


PROIBLEM 47 


If blocks A, B, and C were contiguous in storage, 
the problem would be a simple 3-way merge, and a count of 
the words moved bo block D would control the merge nicely. 
The three blocks could be moved to become contiguous, of 
course, but that would be cowardly and moreover would be 
impractical for longer blocks. The accompanying 
flowchart shows a scheme (by Edward Ryan) for the merge. 


A Merging Problem 


a i a a a 


POPULAR COMPUTING's rating scale for calculators appeared 
in issue No. 10 and was repeated in issue No. 30. The prices 
on pocket calculators are still dropping, and the capabilities 
are still increasing. Since our rating scale is a performance/ 
price ratio, the ratings keep climbing. Consider two widely 
available machines: € 


1. The APF Mark 51 (sold by Sears at $35). This machine 
is an 8-digit floating point calculator for arithmetic operations. 
When it goes into scientific notation, it uses a 5-digit 
mantissa and a e-digit exponent, with a range on the exponent 
of +99. Trig functions, direct and inverse, and natural logs, 
direct and inverse, all operate to 5 significant digits. There 
is one word of accumulating storage that functions either + or -. 
Trig functions can be in degrees or radians, At $35, this 
machine has a rating on our scale of 32.9. 


2. The Texas Instruments! SR-50 is now being discounted 
as low as $75. it is a fully scientific machine, with all 
trig and log functions (except common antilog) and every function 
is calculated to 13 significant digits, with 10 digits in the 
display plus a e-digit exponent. Thus, for all pairs of 
inverse operations (log and exponential; cosine and arccosine; 
and so on) it is difficult to find any number for which such a 
pair of functions does not return to the original value. For 
example, upon entering a number and depressing the reciprocal 
key repeatedly, the original number appears exactly as entered 
every other time. At $75, the SR-50 rates 35.0. 4096 

On an 8-digit machine, the operation (4097/4096) ~ 

--found by squaring 12 times--approximates the natural log base, 
yielding 2.7169. With the extra precision of the SR-50, this € 
trick can be extended to 


(524289/524288 ) 924288 
---squaring 19 times--to yield 2.7182775. L] 
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A subroutine in a program produces as its output 
3-digit numbers, (N), in no special order, and with 
possible strings of duplicates. For example, the 
output on successive calls of the subroutine could be 
the following sequence: 


387 387 123 006 006 527 527 527 527 129 123 123 
387 387 387 387 387 109 109 123 123 123 687 687 
687 687 687 687 687 109 109 109 109 153 153 5e7 
387 387 387 687 123 006 527 109 109 109 109 109 
109 687 109 496 345 456 006 006 006 006 006 006 
006 006 006 109 005 456 567 144 1H4 238 ... ... 


A second subroutine is needed which will examine 
these numbers as they are produced (that is, the second 
subroutine will be called after each new value of N 
is produced) and report each new longer string of 
duplicates. Thus, starting with the sequence of N 
values given above, the printed output from the second 
subroutine should be: 


Bon 2 
5e7 4 
387 5 
687 t 
006 9 and so on. 


A). Flowchart the logic of the second subroutine. 


prosiem 117 


B). Devise a procedure to test that logic. 


[Note: this is Problem 45 (PC13-13) repeated, but in a more explicit 
form. No solution to Problem 45 has been received.] 


Counting String Lengths 


